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in TPD doped Polystyrene 
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Abstract 

"?. 
P3 , A total negative field dependence of hole mobility down to low temperature was ob- 

c/3 I served in N,N'-diphenyl-N,N'-bis(3-methylphenyl)-(l,l'-biphenyl)-4,4'diamine (TPD) 

doped in Polystyrene. The observed field dependence of mobility is explained on the 

basis of low values of energetic and positional disorder present in the sample. The 

low value of disorder is attributed to different morphology of the sample due to 

aggregation/crystallization of TPD. Monte Carlo simulations were also performed 

'^ ■ to understand the infiuence of aggregates on charge transport in disordered medium 

C , with correlated site energies. The simulation supports our experimental observations 

fi I and justification on the basis of low values of disorder parameters. 
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1 Introduction 



Films of organic semiconducting materials are widely used in developing vari- 
ous optoelectronic devices like organic light emitting diodes (OLED), organic 
field effect transistors (OFET), organic solar cells etc [1,2]. Spin cast films 
of molecularly doped polymers (MDP) are used most widely as active layer 
in these devices because of simple fabrication technique. These MDP films 
are mostly amorphous/disordered and have low carrier mobility. To improve 
mobility and transport properties various methods, like annealing, irradia- 
tion etc. are employed to reduce the disorder in the active layer (i.e. change 
film morphology) [3-6] . These methods improve structural order in the organic 
thin film by creating more ordered regions in the otherwise amorphous film. 
Often, depending upon the processing conditions while film casting, the or- 
dered regions are also formed unintentionally due molecular aggregation / 
crystallization of dopants in MDP or aggregation of polymer chains in case 
of conjugated polymers [7,8]. Therefore, in most of the cases where the active 
layers are MDP or pure polymer, films are not purely disordered or amorphous 
rather they are partially ordered. The charge transport in the partially ordered 
films therefore demands deeper understanding of transport mechanism. The 
assumptions made in hopping based charge transport models, for example 
the Gaussian Disorder Model(GDM) [9], developed for a completely isotropic 
and disordered medium probably not sufficient [10]. Moreover these models of 
charge transport do not consider morphology of the active layers in finer de- 
tails. The presence of ordered regions can in-fact reduce the overall energetic 
disorder in the material [5,7] and that enhances the mobility but substantially 
influence the charge transport mechanism [7,11,12]. Field and temperature 
dependence of mobility in such partially ordered samples show wide variation 
with morphology of the sample and is subject of intense research. In general 
the field dependence of the mobility in polymer films follows Poole-Frenkel 
type behavior. However, in certain cases, at low temperature, either very weak 
field dependence or even negative field dependence of mobility has also been 
reported [13-16]. Reports also suggest dependence of mobility on temperature 
as 1/T as well as 1/T^ dependence. In some cases it is hard to distinguish 
between the two power law behavior of the mobility [13,14,17]. 

In Poole-Frenkel type behavior, the mobility increases with increase of electric 
field in a log fiVs.E^^"^ fashion. According to GDM the increase of mobility 
in logfiVs.E^^'^ fashion in disordered molecular solids is due to the tilting of 
density of states by the applied potential that lead to the decrease of energetic 
barrier as seen by the charge carriers in its transit [9]. At higher electric field 
strength the energetic barrier seen by the carrier is negligibly small and this 
results in the saturation of drift velocity of the carrier. Once the drift velocity 
get saturated the mobility decrease with further increase of electric field, i.e. 
negative field dependence of mobility[18]. At intermediate field strength when 



the temperature is low, the field dependence of mobility remains positive. 
When temperature increases carrier gains more thermal energy to over come 
the energetic barrier. Thus the mobility increases and the field dependence of 
mobility weakens. At higher temperature energetic barrier seen by the carrier 
is negligibly small which results in negative field dependence of mobility even 
at intermediate field strengths. Lower the energetic disorder inside the sample 
lower is the temperature at which negative field dependence can be observed 
[7a,ll,13b]. In principle if disorder is very low one can observe negative field 
dependence of mobility at lower temperature. The above reasoning is also 
justified on the basis of GDM. GDM predicts the variation of the mobility as 
given by equation [9], 



fi = HqCxp 
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where /io is a prefactor mobility, a is the measure of energetic disorder (width 
of the Gaussian distribution of site energies), E is the measure of positional 
disorder, the measure of geometrical disorder, k is the Boltzmann constant, T 
is the temperature in Kelvin and C is an empirical constant. From Equation 
(1), the term { (^)^ — S^ !• decides the slope of field dependence of mobility at 

the intermediate field regime, where the mobility shows log nVs.E^'^ depen- 
dence. When temperature increases the term becomes negative. Also, if the 
value of a is low then the above term become negative at lower temperatures. 
This suggests that if the overall energetic disorder of the film is very low or 
if the charge transport occurs through regions of very low energetic disorder 
then the slope of field dependence of mobility can remain negative down to 
lower temperatures. Similarly if (^)^ < S^, the mobility shows a negative 
field dependence. This generally happens only when the positional disorder 
in the sample is remarkably very high. It has also been shown using Monte 
Carlo simulations that high positional disorder can lead to negative field de- 
pendence of mobility at lower electric field strength [9]. In MDPs the negative 
field dependence has generally been observed at very high temperature and 
also when the concentration of the dopant is very low [3] . 

In this paper we show negative field dependence of mobility down to low 
temperatures in MDP. We have investigated the field and temperature de- 
pendence of mobility in films of N,N'-diphenyl-N,N'-bis(3-methylphenyl)-(l,l'- 
biphenyl)-4,4'diamine (TPD) dye. TPD dye is a well known blue emitting laser 
dye and also a hole transporting dye. Films of TPD have been used widely in 
fabricating various organic devices [2,3]. We used TPD doped in Polystyrene 
(PS) at 40:60 proportions by weight (TPD:PS) and measured hole mobil- 
ity using time-of-fiight(TOF) transient photoconductivity technique. A total 
negative field dependence of mobility was observed down to low temperature 
~150K. This was attributed to the low value of disorder in the sample due 



to the presence of aggregation/crystallization of TPD molecules. The study 
highlights the drastic change in morphology of the sample and the resulting 
reduction in the overall disorder of a molecularly doped polymer upon aggre- 
gation/crystallization of dopant. Monte Carlo simulation was performed to 
understand the influence of aggregates on charge transport with correlated 
site energies. Simulation supports our experimental observation as well as the 
justification provided on the basis of low disorder present in the sample. 



2 Experimental and Simulation Details 

2.1 Experimental details 



Solutions of TPD doped Polystyrene at 40:60 proportions by weight were made 
by dissolving the required amount of TPD and PS in Chloroform. Films were 
spin-cast on to a neatly cleaned fluorine doped tin oxide (FTO) coated glass 
substrate. Samples were kept in vacuum for 24 hours at ambient temperature 
to remove the residual solvent. A thin layer of amorphous Selenium (a-Se) and 
an Aluminum top electrode was deposited on to it by thermal evaporation. All 
coatings were done at base pressure of 10~^mbar. Thickness and capacitance 
of the samples were ~4/im and ~10pF respectively. Field and temperature 
dependence of mobility in these samples was determined using conventional 
small signal Time of Flight (TOF) transient photoconductivity technique [3]. 
Samples were mounted on a homemade cryostat to perform temperature de- 
pendent studies. A variable DC potential was applied across the device such 
that no injection occurs from the electrodes to the sample. A 15ns laser pulse 
from second harmonic of Nd: YAG laser (532nm) was used for generating 
a thin sheet of charge in a-Se at a-Se/ TPD:PS interface. Laser intensity is 
adjusted so that total charge generated is less than 0.05 CF, where C is the 
capacitance of the device and V is the voltage applied across the device. The 
time resolved photocurrent is acquired using an oscilloscope as a voltage across 
a load resistance. The transit time, r, is obtained from photocurrent signal 
and the mobility is calculated using /i = L'^/Vt , where L is the thickness 
of the sample. Morphology of the samples was characterized using Scanning 
electron microscope (SEM) images. 



2.2 Details of Monte Carlo Simulation 



Monte Carlo simulations wereperformed to support our experimental obser- 
vations and explore the influence of aggregates on field and temperature de- 
pendence of mobility. The Monte Carlo simulation is based on the commonly 



used algorithm reported by Schonherr et al [19]. A lattice of 70x70x70, along 
X, y and z direction, with lattice constant a = 6A was used for computa- 
tion. Z direction is considered as the direction of the applied field. The size 
of the lattice is judged by taking into account of the available computational 
resources. The site energies of lattice were initially taken randomly from a 
Gaussian distribution of mean ~5.1eV and standard deviation a = 75meV, 
which gives the energetic disorder parameter a=a/kT . The value for a was 
chosen close to the experimental value observed in TPD based MDPs [3]. The 
site energies were made correlated by considering the energy of a site as an 
average of energies of neighboring sites which is defined as follows [20], 

e. = N-'/'J2Kiu,)e, (2) 



Where the variable Ej denotes the uncorrelated energies on the neighboring 
sites, A'^ is the normalization factor that results in required standard deviation 
and K is kernel that provides a degree of correlation among the sites [20] . In 
our simulation, kernel K is considered as unity within a sphere of radius a (a is 
the intersite distance) and zero outside this sphere. Simulation was performed 
on this energetically disordered lattice with the assumption that the hopping 
among the lattice sites was controlled by Miller-Abrahams equation [21] in 
which the jump rate Uij of the charge carrier from the site i to site j is given 
by 



AR,,\ (exp (-2^) exp (±f?) ^ e, > e. 



iy., = Uoexpi-2^a^^ ^ V ..y ^ V - y . (3) 



1 -^ Ej < Ei 



where E is the applied electric field, a is the intersite distance, k is the Boltz- 
mann constant, T is the temperature in Kelvin, Ai?jj = Ri — Rj is the distance 
between sites i and j and 27a is the wave function overlap parameter which 
controls the electronic exchange interaction between sites. Throughout the 
simulation we assume 27a = 10[9,19]. 

Film morphology was varied by incorporating cuboids of so called ordered re- 
gions (representing molecular aggregates/microcrystallites in MDPs) of vary- 
ing size that are placed randomly inside the otherwise disordered host lattice. 
Sizes of ordered regions were limited to a maximum size of 25x25x40 sites along 
X, y and z directions. The energetic disorder inside the ordered region was kept 
low compared to the lattice. This is justified because the aggregates are more 
ordered regions and hence to simulate the charge transport the cuboids must 
be of low energetic disorder compared to host lattice. The site energies inside 
the ordered regions were also taken randomly from another Gaussian distribu- 
tion of standard deviation ~15meV (we chose 5 times less compared to host 
lattice) . Earlier reports have even suggested a ten fold reduction of energetic 



disorder inside the aggregates [22]. Further details of simulation with varying 
film morphology are provided in Ref [1 1] . The site energy of the ordered region 
was also correlated as explained above. The mean energies of ordered regions 
were chosen such that their difference from the mean energy of host lattice is 
in the order of kT. This is justified by the fact that the aggregation of dopants 
can also lead to change in the energy gap (shift in HOMO, LUMO levels) and 
hence the mean energy of Gaussian distribution. Simulations were performed 
by varying the concentration of such ordered regions (varying the percentage 
of volume of lattice occupied by ordered region), temperature and electric field 
so as to simulate the field and temperature dependence of mobility. 



3 Results and Discussions 



3. 1 Experimental Results 



Fig.l shows the typical time of flight transient signal obtained for TPD:PS 
(40:60 wt%) at 290K. Transient showed some plateau suggesting non-dispersive 
transport. We observe dispersive transport at low electric field strength and 
low temperature. The transit time was always determined from the double 
logarithmic plot as shown in the inset of Fig. 1. Fig. 2 shows the field depen- 
dence of mobility parametric with temperature. Negative field dependence of 
mobility was observed down to low temperature and through out the range of 
electric field studied. The mobility becomes almost field independent at 150K. 
There is no remarkable increase in mobility above room temperature compared 
to lower temperatures. Data at low electric field strength, at low temperature, 
were not recorded due to very low magnitude of photocurrent transient signal. 
At high temperatures, in the low field regime, a slight increase in mobility was 
observed with decrease of electric field. Temperature dependence of zero field 
mobility and the slope {(3) of intermediate field region of \og jji{E = 0)Vs.E^^'^ 
follow T~^ as predicted by GDM (Figure 2(b) and 2(c) respectively) [9]. Based 
on GDM formalism we estimated the disorder parameters. The obtained values 
of a, S and fio are 0.0287eV, 2.277 and 2.8471x10^2 cmVVs respectively. The 
value of energetic disorder and positional disorder in our sample is low com- 
pared to the reported values for similar and other molecularly doped polymer 
systems [3,9]. Hence the observed negative field dependence of mobility down 
to low temperature (~150K) in TPD:PS can be attributed to the presence of 
very low value of energetic disorder in the sample studied. As explained above 
(see introduction), low of value of energetic disorder can lead to negative field 
dependence of mobility down to low temperature. When the disorder is very 
low then the barrier offered to carrier is also very low. Thus drift velocity of 
the carrier saturates at lower electric field strength which lead to the decrease 
of mobility with increases of field strength [9,16]. 



The morphology of TPD:PS in our study was investigated using SEM images. 
Fig. 3 shows the SEM images of TPD:PS in our study. SEM images showed 
that TPD has undergone aggregation. Aggregation of TPD resuhed in the 
formation of crystalhne regions of few micron sizes and also chains of such 
microcrystals that are spread all over the sample. Earlier reports of crystal- 
lization/aggregation of TPD in TPD doped Polystyrene system, even at low 
concentration and without annealing [23], also supports our observation. Ear- 
lier reports of mobility measurement, who apparently report only positive field 
dependence, in TPD doped polymers assert that their study was performed 
in totally amorphous sample with no signs of aggregation of dopants at least 
during the time of experiment [2,3,24]. Presence of more ordered regions can 
drastically change the entire morphology of the sample and can lead to reduc- 
tion in the effective energetic and positional disorder. This is consistent with 
our observation of low energetic and positional disorder in TPD:PS. Due to 
aggregation the charge transport therefore occurs through a combination of 
ordered and less ordered regions. In such cases the charge transport is highly 
influenced by the packing and orientation of ordered regions. If the charge 
transport occurs mostly through these ordered regions, regions of very low 
disorder, then a negative field dependence of mobility down to low tempera- 
ture can be expected. So the presence of these microcrystallites can drastically 
change the morphology with lower energetic and positional disorder and the 
behavior of charge transport in these samples. 



3.2 Simulation Results 



In order to justify and support the above explanation a Monte Carlo simula- 
tion of charge transport in disordered lattice was also performed. Aim of the 
simulation was to understand the influence of aggregates on charge transport, 
in particular the negative field dependence. Simulated field and temperature 
dependence of mobility for a pure host lattice having DOS with standard 
deviation ~75meV (without considering any ordered regions) has been dis- 
cussed in detail in our earlier report [11]. Simulation results were as predicted 
by GDM. To study the influence of embedded micro crystals/aggregation of 
dopants on charge transport, the simulation was performed after incorporating 
ordered regions inside the host lattice (as explained in simulation procedure 
above). Fig. 4 shows the field dependence of mobility, at ~248K, parametric 
with the concentration of ordered regions having DOS with standard devia- 
tion ~15meV and mean energy lower by ~ kT compared to mean energy of 
host lattice. Magnitude of mobility at all regimes of electric field increases 
with increase of concentration of ordered regions concomitant with decrease 
of slope at the intermediate field regime. The saturation of mobility and de- 
crease of mobility with further increase of electric field, at high field regime, 
was observed when the concentration of ordered region was higher than 60%. 



With increase in concentration of ordered region the saturation of mobihty 
occurs at lower electric field strength. For the cases when the concentration 
of ordered region is less that 60% the mobility was not completely saturated 
even for the maximum strength of electric field used for simulation. When 
the concentration of the ordered regions inside the host lattice is very high 
(~98%) the field dependence of mobility become totally negative even at this 
low temperature (~248K). The mobility decreases with increase of electric 
field right from low electric field strengths used for simulation. 

The observed features in the field dependence of the mobility, after incorpo- 
rating the ordered regions inside the host lattice, can be explained on the basis 
of decrease of effective energetic disorder in the host lattice when ordered re- 
gions are embedded in it. If the energetic disorder is small then the energetic 
barrier seen by the carriers due to disorder will be small. This results in higher 
mobility but weaker field dependence. This explains the increase in magnitude 
of mobility and decrease in the slope of semi log fiVs.E^^"^ curve at intermedi- 
ate field regime with increase in concentration of ordered regions in the host 
lattice. The low value of energetic disorder also can lead to saturation of drift 
velocity /mobility at lower electric field strength. When the concentration of 
ordered regions inside the host lattice is very high, the overall energetic disor- 
der will be very low. At such high concentration the carrier will travel mostly 
through ordered region only. Since the energetic disorder inside the ordered re- 
gion is small (only ~15meV in our study) the mobility saturates at low electric 
field strength and even at low temperatures. In such a case one can observe 
total negative field dependence of mobility down to low temperatures. More- 
over for a very low value of energetic disorder the term I {-^Y ~ ^^ f remains 
negative down to low temperature. It is also possible to have total negative 
field dependence down to low temperature when the charge transport occurs 
totally through ordered regions, i.e. when the whole host lattice is completely 
occupied with ordered region (Fig. 4, 100%). Fig. 5 shows the field dependence 
of mobility parametric with temperature when host lattice is completely occu- 
pied with ordered region. A total negative field dependence of mobility down 
to low temperature is observed but the mobility in this case decreases with 
increase of temperature. Higher magnitude of mobility was observed for lower 
temperature (~150K) and low magnitude of mobility at higher temperature 
(~390K). This kind of behavior is possible when the overall energetic disorder 
is very low. When energetic disorder is very low and the thermal energy is 
comparable to energetic disorder then thermal energy dominates which forces 
the carriers to move in longer paths. Effectively the transit time increases and 
mobility decreases with increase of temperature [16]. 

From the simulation it is inferred that charge transport in TPD:PS occurs 
mostly through aggregates but not completely through aggregates. Hence the 
effective disorder seen by the carrier is low which results in total negative 
field dependence of mobility. Experimentally we observed that mobility de- 



creases with decreases of temperature which rules out the possibihty that 
charge transport is occurring completely through aggregates. The influence of 
small disordered regions present in entire trajectory of charge transport (major 
part of the trajectory is occupied by ordered region) for each carrier become 
prominent only at low temperature. When the thermal energy of carrier is 
low, the carrier will see the effect of disorder and that results in the positive 
field dependence of mobility at low temperature. Fig. 6(a) shows the field 
dependence of mobility parametric with temperature when 98% of the host 
lattice is occupied with ordered regions having DOS with standard deviation 
~15meV and mean energy lower by ~ kT compared to mean energy of host 
lattice. Mobility decreases with decrease of temperature and positive field de- 
pendence of mobility was observed at lower temperature (200K). Negative field 
dependence of mobility was observed down to 248K. At temperature around 
300K, there is no variation of mobility with increase of temperature for all 
electric field strength used for simulation. This observation is similar to what 
observed experimentally in TPD:PS at higher temperatures. The temperature 
dependence of zero field mobility and slope of intermediate field region of log 
versus E^^'^ plot followed 1/T^ dependence as predicted by GDM (Data not 
shown) . 



4 Conclusion 



Time of flight mobility measurement was carried out in TPD:PS. A total nega- 
tive field dependence of mobility down to low temperature was observed. This 
was explained on the basis of low value of energetic disorder in TPD:PS. The 
low value of energetic disorder was attributed to aggregation of TPD that re- 
sults in the formation of crystalline regions of few micron size and chaining of 
such crystals. The presence of aggregates drastically changes the morphology 
of the sample and reduces the overall disorder in the sample. Monte Carlo 
simulation study clearly showed the influence of aggregates on charge trans- 
port which supports our experimental observation. The simulation suggests if 
charge transport occurs mostly through aggregates (crystalline regions of very 
low energetic disorder) the mobility saturates at low electric field strengths and 
it can lead to a total negative field dependence of mobility down to low temper- 
ature. From the simulation studies it is inferred that the charge transport in 
TPD:PS in the present study are not occurring totally through aggregates but 
in combination of aggregates and disordered medium. Thus our experimental 
and simulation studies showed the influence of morphology of the sample on 
field and temperature dependence of mobility. The study also highlights the 
need to consider film morphology in various models used for analyzing charge 
transport in polymer films. 
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Figure Captions 



Figure 1. Typical time of flight transient signal for TPD:PS at 4xlO^V/cm, 
290K. Inset show the double logarithmic plot used for determining the transit 
time. 

Figure 2. Experimental data for TPD:PS(40:60) (a) Electric field dependence 
of mobility semi log^Vs-E^^'^ parametric with temperature, (b) semi log iJ,{E = 
0) Vs. T~^ (c)/9 Vs. {a /kTY plot, where (3 is the slope of high field region of 
sem,i\ogiiVs.E^/'^ plot. 

Figure 3. SEM image of TPD:PS showing crystallization/aggregation of TPD 
molecules. 

Figure 4. Simulated field dependence of mobility for host lattice (o'=75meV), 
at 248K, embedded with ordered regions of various concentrations having DOS 
of standard deviation~15meV. 

Figure 5. Field dependence of mobility of host lattice with ~15meV paramet- 
ric with temperature. Case is similar when host lattice is completely occupied 
with ordered regions. 

Figure 6. Field dependence of mobility of host lattice, parametric with tem- 
perature, with embedded ordered regions having DOS with standard deviation 
~15meV. Concentration of embedded ordered region inside the host lattice is 
~98%. Straight line show linear fit to intermediate field region. 
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Figure 2(a) 
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Figure 2(b) 
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Figure 2(c) 
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Figure 3 
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Figure 4 
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Figure 6 
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